General construction of noiseless networks detecting entanglement 

with help of linear maps 
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We present the general scheme for construction of noiseless networks detecting entanglement 
with the help of linear, hermiticity-preserving maps. We show how to apply the method to detect 
entanglement of unknown state without its prior reconstruction. In particular, we prove there always 
exists noiseless network detecting entanglement with the help of positive, but not completely positive 
maps. Then the generalization of the method to the case of entanglement detection with arbitrary, 
not necessarily hermiticity-preserving, linear contractions on product states is presented. 
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I. INTRODUCTION 

It has been known that entanglement can be detected 
with help of special class of maps called positive maps 
[0> S HI- In particular there is an important criterion 
l] saying that g acting on a given product Hilbcrt space 
TLa ® T~(-b is separable if and only if for all positive (but 
not completely positive) maps A : B(Hb) — * B(Ha) 0] 
the following operator 



X A (g) = [I ® A]( e ) 



(1) 



has all non-negative eigenvalues which usually is written 
as 



[I® A](g) > o. 



(2) 



Here by I we denote the identity map acting on B(Ha)- 
Since any positivity-preserving map is also hermiticity- 
preserving, it makes sense to speak about eigenvalues of 
X\(g). However, it should be emphasized that there are 
many As (and equivalently the corresponding criteria) 
and to characterize them is a hard and still unsolved 
problem (see, e.g., Ref. [f| and references therein). 

For a long time the above criterion has been treated as 
purely mathematical. One used to take matrix g (ob- 
tained in some prior state estimation procedure) and 
then put it into the formula @. Then its spectrum was 
calculated and the conclusion was drawn. However it can 
be seen that for, say states acting on H.a®'H.b ~ C d £g>C d 
and maps A : B(C d ) — ► B(C d ), the spectrum of the opera- 
tor X A (g) consists of n spoc = d 2 elements, while full prior 
estimation of such states corresponds to n cst = c? 4 — 1 pa- 
rameters. 

The question was raised @ as to whether one can per- 
form the test ^ physically without necessity of prior 
tomography of the state g despite the fact that the map 
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7® A is not physically realizable. The corresponding an- 
swer was @ that one can use the notion of structural 

physical approximation I ® A (SPA) of un-physical map 
i" ® A which is physically realizable already, but at the 
same time the spectrum of the state 



X A (g) = {I ® A](g) 



(3) 



is just an affine transformation of that of the (unphysical) 
operator Xa(q). The spectrum of Xa(q) can be measured 
with help of the spectrum estimator [7|, which requires 
estimation of only d 2 parameters which (because of affin- 
ity) are in one to one correspondence with the needed 
spectrum of 0. Note that for 2 <8> 2 systems (the com- 
posite system of two qubits), similar approaches lead to 
the method of detection of entanglement measures (con- 
currence @ and entanglement of formation Q) without 
the state reconstruction [10 ]. 

The disadvantage of the above method is [ll[ that re- 
alization of SPA requires addition the noise to the system 
(we have to put some controlled ancillas, couple the sys- 
tem, and then trace them out). In Ref. [ll| the question 
was raised about the existence of noiseless quantum net- 
works, i.e., those of which the only input data are: (i) 
unknown quantum information represented by g® m (ii) 
the controlled measured qubit which reproduces us the 
spectrum moments (see Ref. 0]). It was shown that for 
at least one positive map (transposition) T the noiseless 
network exists [ill ]. Such networks for two-qubit con- 
currence and three-qubit tangle have also been designed 

In the present paper we ask a general question: do 
noiseless networks work only for special maps (functions) 
or do they exist for any positive map test? In the case 
of a positive answer to the latter: is it possible to de- 
sign a general method for constructing them? Can it be 
adopted to any criteria other than the one defined in |(5J? 

For this purpose we first show how to measure a spec- 
trum of the matrix 0{g), where 6 : B{C m ) -► B(C m ) 
is an arbitrary linear, hermiticity-preserving map and 
g is a given density operator acting on C m , with the 
help of only m parameters estimated instead of m 2 — 1. 
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For bipartite g where m — d 2 this gives d 2 instead of 
d — 1. This approach is consistent with previous results 
, [lH where arbitrary polynomials of elements of a 
given state g have been considered. In these works it was 
shown out that any at most k-th degree polynomial of a 
density matrix g can be measured with help of two col- 
lective observables on k copies of g. In fact one can treat 
the moments of Q(g) which we analyze below as polyno- 
mials belonging to such a class. We derive the explicit 
form of observables for the sake of possible future ap- 
plication. Moreover, approach presented in the present 
paper allows for quite natural identification of observ- 
able that detects an arbitrary polynomial of the state 
g subjected to some transformation 9. Then we pro- 
vide an immediate application in entanglement detection 
showing that for suitable O the scheme constitutes just 
a right method for detecting entanglement without prior 
state reconstruction with the help of cither positive map 
criteria ^ or linear contraction methods discussed later. 



II. GENERAL SCHEME FOR CONSTRUCTION 
OF NOISELESS NETWORK DETECTING 
SPECTRUM OF @(g) 

A. Construction of an observable 

Since m x m matrix Q(g) is hermitian its spectrum 
may be calculated using only m numbers 



a k = Tr[9(p)f = 



A 



(fc = l,...,m), (4) 



but permutes states |e») in the reversed order. Therefore 
all the numbers a k may be expressed as 



1 



a k 



-Tr 



(V (fe) + V< fc >t) 0{g) m ] . (7) 



Let us focus for a while on the map 9. Due to its 
hermiticity-preserving property it may be expressed as 



9( 



E 



(8) 



with r/j € K and Kj being linearly independent m-by-m 
matrices. By the virtue of this fact and some well-known 
properties of the trace, after rather straightforward alge- 
bra we may rewrite Eq. ([7]) as 



1 

oik = — Tr 
2 



(et)®* (y {k) + V {k A g m 



(9) 



where 9^ is a dual map to 9 and is given by 0'(-) = 
Y^iVi^H (')Ki- Here we have applied a map (9 t )® fc on 
the operator V^+V^ instead of applying 9® fc to g® k . 
This apparently purely mathematical trick with the aid 
of the fact that the square brackets in the above contain 
a hermitian operator allows us to express the numbers 
Qfc as a mean value of some observables in the state g® k . 
Indeed, introducing 



( e k) = ^Tr 



(etf k (v^ + v^)} (10) 



we arrive at 



where Xi are eigenvalues of 9(g). We shall show that 
all these spectrum moments can be represented by mean 
values of special observables. To this aim let us consider 
the permutation operator defined by the formula 



V (k) 



ei>|e 2 ; 



\ e k) 



|e fc ) |ei> (g) ... (g) |e fc . 



(5) 



where (k = 1, . . . , m) and | e^} are vectors from C m . One 
can see that V^' is just an identity operator l m acting 
on C m . Combining Eqs. (|4]) and ([5]) we infer that ak 
may be expressed by relation 



o-k = Tr 



{V fe) [9(e)] 0fc } 



(G) 



which is generalization of the formula from Refs. @, 0] 
where 9 was (unlike here) required to be a physical oper- 
ation. At this stage the careful analysis of the right-hand 
side of Eq. ([6]) shows that ak is a polynomial of at most 
fc-th degree in matrix elements of p. This, together with 
the observation of Refs. [H, [TJ, [l5[ allows us already 
to construct a single collective observable that detects 
ak- However, for the sake of possible future applications 
we derive the observable explicitly below. To this aim we 
first notice that ak may be obtained using hermitian con- 
jugation of which again is a permutation operator 



a k = Tr 



/n(*0 „0k 

ft y 



(ii) 



In general, a naive measurement of all mean values 
would require estimation of much more parameters that 
m. But there is a possibility of building a unitary network 
that requires estimation of exactly m parameters using 
the idea that we recall and refine below. 

Finally, let us notice that the above approach general- 
izes measurements of polynomials of elements of g in the 
sense that it shows explicitly how to measure the polyno- 
mials of elements of 9(g). Of course, this is only of rather 
conceptual importance since both issues are mathemati- 
cally equivalent and have the origin in Refs. (l3l. [l5j . 



B. Detecting mean of an observable by 
measurement on a single qubit revised 

Let A be an arbitrary observable (it may be even in- 
finite dimensional) which spectrum lies between finite 



numbers 



and 



and a be a state acting on H. 



In Ref. [lj| it has been pointed out that the mean value 
(A) a = Tr„4er may be estimated in process involving the 
measurement of only one qubit. This fact is in good 
agreement with further proof that single qubits may serve 
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as interfaces connecting quantum devices [lTf. Below we 
recall the mathematical details of the measurement pro- 
posed in Ref. 16]. At the beginning one defines the 
following numbers 



t ' 



t'A 



a {+) = a { - ] 
a A — a A 



a ( X ] = max{0,-a2 in }, 
and observe that the hermitian operators 



and 



V = J[a { A h n +A 



S+) 



(12) 



(13) 



FIG. 1: General scheme of a network for estimating mean 
value of an observable A, with a bounded spectrum, in a given 
state a. Both U-h' and its conjugate U 1 ^ standing before Ua 
can obviously be removed as they give rise to identity, last 
unitary on the bottom wire can be removed as it does not 
impact measurement statistics on the top qubit. However, 
they have been put to simplify subsequent network structure. 



Vi 



VjVo 



(14) 



satisfy J2i=o = ^-H EH and as such define a general- 
ized quantum measurement which can easily be extended 
to a unitary evolution (see Appendix A of Ref. [13] for a 
detailed description). Consider a partial isometry on the 
Hilbert space C 2 ® H defined by the formula 



l 

^ = 5>>(0|®V- 

1=0 



V 
V x 



(15) 



The first Hilbert space C 2 represents the qubit which 
shall be measured in order to estimate the mean value 
(A) a - The partial isometry can always be extended to 
unitary Ua such that if it acts on |0)(0| ®u then the final 
measurement of observable cr 2 [19] on the first (qubit) 
system gives probabilities "spin- up" (of finding it in the 
state |0)) and "spin-down" (of finding in state |1)), re- 
spectively of the form 



Po = Tr (VJV 0Q ) , 



= Tr(v?V ig ) = l-p . (16) 



Pi 



One of the possible extensions of U A to the unitary on 
C 2 ® TC is the following 



V -Vx 
Vx V 



Vn - icr„ 



Vx. 



(17) 



where Vo is a projector onto state |0), i.e., Vo — |0)(0|, 
then we have an easy formula for the mean value of the 
initial observable A: 



(A), 



= a^po 



J-) _ a (+)M±l 

a A — a A 2 



-,(-) 



(20) 



A general scheme of a network estimating the mean 
value (|2T)]) is provided in Fig. [1] We put an additional 
unitary operation on the bottom wire after unitary Ua 
(which does not change the statistics of the measurement 
on control qubit) and divided identity operator into two 
unitaries acting on that wire which explicitly shows how 
simplification introduced in Eq. (I18j) works in practice. 

Now one may ask if the mean value (A) a belongs to 
some fixed interval, i.e., 



cx < (A)a < c 2 , 



(21) 



where ci and ci are real numbers belonging to the spec- 



trum of A, i.e., [a™ : 



x ] (e.g. if A is an entanglement 



witness and we want to check the entanglement of a state 
a then we can put c\ = and ci — a^J ax , and condition 
(f2"Tj) reduces to (A) a > 0). Then one easily infers that 
the condition (f2"Tj) rewritten for visibility is 



,(-) 



,(+) 



1 < v A < 2 



C2 



,(-) 



J A 



1. 



(22) 



The unitarity of U A follows from the fact that operators 
Vo and V± commute. Due to the practical reasons in- 
stead of unitary operation representing POVM {Vo, Vx} 
we shall consider 

U dct (A, \J' H ) = (l a <8> XJ' n ) U A (1 2 ® , (18) 

where I2 is an identity operator on the one-qubit Hilbert 
space C 2 and is an arbitrary unitary operation that 
acts on H and simplifies the decomposition of U A into 
elementary gates. Now if we define a mean value of mea- 
surement of a z on the first qubit after action of the net- 
work (which sometimes may be called visibility): 



va 



Tr 



{<Jz ® In) (I2 ® U H ) U A V 



■aU A 



(19) 



Having the general network estimating va, one needs 
to decompose an isometry Ua onto elementary gates. 
One of possible ways to achieve this goal is, as we shall 
see below, to diagonalize the operator Vq. Hence we may 
see Eq. QUI) to be 



choose U H 



(23) 



with |<^k) being normalized eigenvectors of Vo indexed 
by a binary number with length 2 k . Since Vo and Vx 
commutes, this operation diagonalizes Vx as well. By 
virtue of these facts, Eq. (fl"g|) reduces to 



U det (A,U w ) = ^C/ k ® |k)(k| 



(24) 
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] Vp*-i — v^-i — \Uy*- 3 \ \W\ — \Ui 



FIG. 2: Noiseless network for estimating moments of O(p) 
with q being a bipartite mixed state, i.e., density matrix act- 
ing on C 2 ® C 2 . 



with unitaries (as previously indexed by a binary num- 
ber) 



Akl2 - iyl - Ako-y 



(25) 



where Ak are eigenvalues of Vo- So in fact we have a 
combination of operations on the first qubit controlled by 
2 k wires. All this combined gives us the network shown 
in the Fig. U 

Now we are in the position to combine all the elements 
presented so far an show how, if put together, they pro- 
vide the general scheme for constructing noiseless net- 
work for spectrum of Q(g) for a given quantum state g. 
For the sake of clarity below we itemize all steps neces- 
sary to obtain the spectrum of 0(g): 



(i) Take all observables 0<#> (k 
Eq. (HO]). 



1. 



, to) defined by 



(ii) Construct unitary operations Uq(u) according to the 

the given prescription. Consider the unitary op- 
eration U det (A, Uft) (JJ't-i arbitrary). Find decom- 
position of the operation into elementary quantum 
gates and minimize the number of gates in the de- 
composition with respect to U^. Build the (opti- 
mal) network found in this way. 

(iii) Act with the network on initial state Vq <g> g® k . 



(iv) Measure the ,, visibilities" v 
cording to (fT9|) . 



(k = 1, . . . , to) ac- 



(v) Using Eq. (j20|) calculate the values of otk (k 
1, . . . , to) representing the moments of ®{g). 



C. Detecting entanglement with networks: 
example 

The first obvious application of the presented scheme 
is entanglement detection via positive but not com- 
pletely positive maps. In fact for any bipartite state 
g € B(Ha <8> Hb) we only need to substitute with 
tj\ ® Ab with Ag being some positive map. Then appli- 
cation of the above scheme immediately reproduces all 



the results of the schemes from Ref. [6j but without ad- 
ditional noise added (presence of which required more 
precision in measurement of visibility) . 

As an illustrative example consider Ab = T, i.e., 
is partial transposition on the second subsystem (usually 
denoted by Tb or by T), in 2 <g) 2 systems. Due to to 
the fact that partial transposition is trace-preserving we 
need only three numbers a^, (fc = 2,3,4) measurable via 
observables 



and 



O 



;3,4) 



o 



(3,4) 



(2) 



(2) 



v 



(3,4)f 



V, 



(2) 



(3,4)t 



V, 



(3,4; 



(26) 



(27) 



where subscripts mean that we exchange first and second 
subsystems respectively. The hermitian conjugation in 
the above may be replaced by transposition since the per- 
mutation operators have real entries. For simplicity we 
show only the network measuring second moment of g Tu . 
General scheme from Fig. [2j reduces then to the scheme 
from Fig. [3] Note that the network can also be regarded 



ict„ H -vjy - -i(j y H -iff., 




FIG. 3: Network estimating the second moment of partially 
transposed two-qubit density matrix g. Uw is decomposed 
to single qubit gates; here U = (l/v / 2)(l2 + icy)- 

as a one measuring purity of a state as Tr(g Tfi ) 2 = Tug 2 . 
Note that the this network is not optimal since an alter- 
native network 0] measuring Trg 2 requires two controlled 
swaps. 



III. EXTENSION TO LINEAR CONTRACTIONS 
CRITERIA 

The above approach may be generalized to the so- 
called linear contractions criteria. To see this let us recall 
that the powerful criterion called computable cross norm 
(CCN) or matrix realignment criterion has recently been 
introduced [U [22| ■ This criterion is easy to apply (in- 
volves simple permutation of matrix elements) and has 
been shown [2l|, [53] to be independent on a positive par- 
tial transposition (PPT) test @]. It has been further gen- 
eralized to the linear contractions criterion [23] which 
we shall recall below. If by qa ( (i = 1, . . . , n) we de- 
note density matrices acting on Hilbert spaces Ti-Ai an d 
by Ti certain Hilbert space, then for some linear map 
K : B[H Al ® . . . <8> H A J -> B(H) we have the following 
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Theorem [23| . If some TZ satisfies 

\\K(QA 1 ®eA 2 ®---®QA n )\\ Tl <l, (28) 

then for any separable state QA\A 3 ...A n € B{HAi ® • • • ® 
Ha„) one has 

\\K(Q A ,A 2 ...A n )\\Tr<l. (29) 

The maps 7?. satisfying (|28|) are linear contractions on 
product states and hereafter they shall be called, in brief, 
linear contractions. In particular, the separability condi- 
tion (j29f comprises the generalization of the realignment 
test to permutation criteria [23!, 01 (see also Ref . [25[ ) . 

The noisy network for entanglement detection with the 
help of the latter have been proposed in Ref. .26] . Here 
we improve this result in two ways, namely, by taking 
into account all maps 1Z of type (f28|) (not only permuta- 
tion maps) and introducing the corresponding noiseless 
networks instead of noisy ones. For the se p urposes we 
need to generalize the lemma from Ref. [26] formulated 
previously only for real maps S : B(TL) — > BiTL). Wc 
represent action of S on any g S B(Tt) as 

S{g) = S^kiTrigP^Pki, (30) 

ij,kl 

where in Dirac notation P xy = \x){y\. Let us define com- 
plex conjugate of the map S via complex conjugation of 
its elements, i.e., 

S*(g) = J2s* jM Tr(gP ij )P kl , (31) 

ij.kl 

where asterisk stands for the complex conjugation. The 
we have the following lemma which is easy to proof by 
inspection: 

Lemma. Let S be an arbitrary linear map on 33(Ti) . 
Then the map S' = [To S* o T] satisfies S'(g) — [S(g)]^ . 

Now let us come to the initial problem of this section. 
Suppose then we have TZ satisfying Eq. ([28]) and a given 
physical source producing copies of a system in state g 
for which we would like to check Eq. (|2"9"]1 . Let us observe 
that 

l|fcG?)lk=EV^ ( 32 ) 

i 

where {7^} are eigenvalues of the operator X-jz(g) = 
TZ(g)TZ(g)^ . Below we show how to find the spectrum 
We need to apply our previous scheme from Sec. |TT] 
to the special case. Let us define the map L n — TZ® TZ', 
where TZ is our linear contraction and TZ' is defined ac- 
cording to the prescription given in the Lemma above, 



i.e., TZ' = [T o TZ* o T]. Let us also put g' = g® 2 and ap- 
ply the scheme presented above to detect the spectrum 
of L-jz(g'). It is easy to see that the moments detected in 
that way are 

Tr[L n (g')] k = Tr [TZ{g)TZ{g)^} * = £ <y* (33) 

i 

From the moments one easily reconstructs {7^} and may 
check the violation of Eq. (|29|) . 

IV. SUMMARY 

We have shown how to detect the spectrum of the oper- 
ator 0(g) for arbitrary linear hermiticity-preserving map 
given the source producing copies of the system in state 
g. The network involved in the measurement is noiseless 
in the sense of [llj and the measurement is required only 
on the controlled qubit. Further we have shown how to 
apply the method to provide general noiseless network 
scheme of detection detecting entanglement with the help 
of criteria belonging to one of two classes, namely, those 
involving positive maps and applying linear contractions 
on product states. 

The structure of the proposed networks is not optimal 
and needs further investigations. Here however we have 
been interested in quite a fundamental question which is 
interesting by itself: Is it possible to get noiseless net- 
works schemes for any criterion from one of the above 
classes? Up to now their existence was known only for 
special case of positive partial transpose (cf. p~2[). Here 
we have provided a positive answer to the question. 

Finally, let us note that the above approach can be 
viewed as an application of collective observables [see Eq 
(fTTj) ]. The general paradigm initiated in Refs. \U1 |27I | 
has been recently fruitfully applied in the context of gen- 
eral concurrence estimates [28], |29| which has been even 
preliminarily experimentally illustrated. Moreover, re- 
cently the universal collective observable detec ting any 
two-qubit entanglement has been constructed [301 ] . It 
seems that the present approach needs further analysis 
from the point of view of collective observables including 
especially collective entanglement witness (see [27I l29l|). 
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